arXiv:1501.07594vl [cs.NI] 29 Jan 2015 


Analytical Model for IEEE 802.15.4 Multi-Hop 
Networks with Improved Handling of 
Acknowledgements and Retransmissions 


Florian Meier 
Institute of Telematics 
Hamburg University of Technology 
Hamburg, Germany 
florian. meier @ tuhh.de 


Volker Turau 
Institute of Telematics 
Hamburg University of Technology 
Hamburg, Germany 
turau@tuhh.de 


Abstract —The IEEE 802.15.4 standard allows for the deploy¬ 
ment of cost-effective and energy-efficient multi-hop networks. 
This document features an in-depth presentation of an analytical 
model for assessing the performance of such networks. It con¬ 
siders a generic, static topology with Poisson distributed data- 
collection as well as data-dissemination traffic. The unslotted 
CSMA/CA MAC layer of IEEE 802.15.4 is closely modeled as 
well as an enhanced model of the neighborhood allows for 
consideration of collisions of packets including interferences 
with acknowledgements. The hidden node problem is taken into 
account as well as a formerly disregarded effect of repeated 
collisions of retransmissions. The model has been shown to 
be suitable to estimate the capacity of large-scale multi-hop 
networks. 
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I. Introduction 


Analytical models have proven to be very suitable and 
effective for modeling large scale wireless mesh networks [1], 
This document is meant to present the model used in this paper 
in detail to allow for better understanding of the underlying 
mathematics and better reproduction of the results. As further 
reference, the reader might also be interested in the work by 


Di Marco et al. [2] that served as a starting point to build this 
model. The following contributions were newly established in 
the model presented below: 

• Downstream traffic from the central gateway to the nodes. 

• Collisions of acknowledgements with packets or other 
acknowledgements by enhancing the conflict graph. 

• Revised computation of collision probabilities after fail¬ 
ures, taking into account simultaneous retransmissions. 

• Minor enhancements to allow for faster computation. 

The source code of the implementation for reproduction 
and further use is published at GitHub [3], 

Fig. 1 shows the different involved modules of the analyti¬ 
cal model and how they are linked. The topology as presented 
in Sect. II is given as an input to the model and describes how 
the nodes are linked and how they route the traffic. These 
information are used for calculating how the actual traffic 
is distributed within the network (Sect. IV) as well as for 
determining which transmissions might lead to busy channels 
and packet collisions in the neighborhood model (Sect. VII). 



Fig. 1: Relationships between the different parts of the model. 


The traffic distribution yields the probability qi that a node 
has a pending transmission to send over link l. Though, it does 
not directly send the packet, but uses a CSMA/CA technique 
with retransmissions to increase the success probability. This is 
modeled in Sect. V and yields the probability t; that it actually 
senses the channel to start the transmission. For this, it also 
depends on information from the neighborhood model, such 
as the probability that the channel is sensed free ai and the 
probability that no acknowledgement arrives and a 

retransmission should be started. 
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The information from the neighborhood model and the 
topology is also used to calculate the final reliability of a 
link Hi, that is the probability that a packet will eventually be 
transmitted over the link and not dropped because of a busy 
channel or repeated collisions. Ri is again used to calculate the 
traffic distribution. This closes a circle, so it is evident that the 
equations of the outer modules are interlinked. They built up a 
non-linear equation system that can be computed numerically. 
For this task, the implementation utilizes the PETSc framework 
[4], [5], The topology, and thereby also the analog model, are 
not in the loop, so they can be calculated offline before the 
actual computation. 


II. Topology 

The model considers a multi-hop network as depicted in 
Fig. 2, with N nodes including the gateway r and N — 1 
clients. The dashed and continuous lines in Fig. 2 depict the 
predicate D{v, w) that a signal sent out by v is strong enough 
to disturb an ongoing reception at w. This boolean predicate is 
calculated by an analog model together with the estimated bit 
error rate BER( t , ^ of a data transmission between v and w. 



Fig. 2: A multi-hop network with a central gateway. 


An exemplary analog model for IEEE 802.15.4 networks 
that is used in the implementation is given in the following. 
It is presented in the IEEE 802.15.4 standard Annex E [6] 
and is based on a breakpoint log-distance path loss model 
determined by measurements. The received signal power in 
dB with distance d V)W and transmission power P tx is given by 

p , \_ P _ / 58.5 + 331og 10 (<4 iW /8rn) d VfW >8m 

r rx[ v,w)-^ tx | 40 2 + 201ogio(d?;iw) d v w <Sm' 

( 1 ) 

It also includes the computation of the bit error rate dependent 
on the noise power Pjy 

SNR ( „ iU)) = IOto'CAA^I-Av^ ( 2 ) 

BHW, = ^4 E (3) 

k—2 ' ' 

and thereby the packet error rate for transmissions of b bytes 

PER 6j(UiW) = 1 ^ (1 - BER;) 8 ' 6 . (4) 


The predicate D(v,w) can be calculated from this, assuming 
there is a minimum interference power Pdist that is strong 
enough to disturb an ongoing reception at w 

D(v,w) ■■= P rx (v,w) > Pdi S t ■ (5) 

In a further preprocessing step, an optimal static routing 
tree T is computed by using Dijkstra’s algorithm. The weight 

— l°g(l — BER (t , )Uj) ) + 10 -3 (6) 

was used in our experiments to minimize the bit error along the 
path, but still minimize the hop count if the BER is negligible. 
Though, any procedure that yields a suitable routing tree T can 
be used. The tree is described by the predicate T(p, c), that is 
true if and only if p is a parent of c in T , yielding the set of 
active links 


L 


up 


L — L„p U L(lown • 


{(c,p) \T{p,c)}, L down := {( p,c ) I T(p,c)}. 


( 7 ) 

( 8 ) 


III. Traffic Generation 


A Poisson traffic generation model is applied, so the 
intervals between two packet generations are exponentially 
distributed. The mean interval for traffic from a node to the 
gateway (upstream) is denoted as / up , while for each client, the 
gateway generates packets with mean interval /down- Analog 
to the IEEE 802.15.4 standard [6], all times are defined as 
multiples of 

Sb ■= aUnitBackoffPeriod • Symbol duration = 20-16 ps. 

( 9 ) 


Taking into account 4 bits per symbol for the O-QPSK PHY, 
a transmission of a packet of length B bytes takes 


8 • B ■ 16 ps B 
A- S b = To 


( 10 ) 


time units. An acknowledgment has a fixed length of B ACK = 
11 bytes. Therefore it takes 

Lack — ^ ( H ) 

time units. Together with the specified supplementary values 
IFS = 40 and t ac k = 12, an acknowledged transmission takes 

IFS + t ac h 

L s = L p + L ACK + -+ L ACK + 2.6 (12) 

time units. The number of symbols t m ^ ac k to wait for a lost ac¬ 
knowledgement ( macAckWaitDuration ) is calculated according 
to the standard 


tm,ack = aUnitBackoffPeriod 

+ aTurnaroundTime + phySHRDuration 
+ [6 • phySymbolsPerOctet ] 

= 20 + 12 + 10 + 6 • 2 = 54. 


( 13 ) 


This gives a total time for an unsuccessful transmission of 

L c = L p + =L p + 2.7. (14) 

time units. In relation to Sb, the packet generation rate on each 
client in upstream direction is 

_ s b 

9u p — T 

-*■ UD 


2 


( 15 ) 



and since the gateway r generates packets for each client its 
packet generation rate is 


//down 


(N - 1) • S b 

^down 


(16) 


IV. Traffic Distribution 

The distribution of the traffic is directly calculated from 
the tree. A/ denotes the packet sending rate of the sender v 
of a link l = (v,w) £ L, including generated and forwarded 
traffic. At the receiver side /j ( is the effective packet rate to 
be forwarded by w. Furthermore, 7 „denotes the number of 
proper descendants of a nodes, that is the number of nodes 
whose traffic is routed through this node (Fig. 3). 




Fig. 3: An example for the set of active links and the values 
of 7 „in the nodes. 


Fig. 4 depicts the possible inflows and drains of traffic. 
For upstream, A ( v , w ) is the sum of the traffic of all proper 
descendants and the traffic generated in this node. Some 
packets will be dropped with a probability of R( v , w ), because 
it was not possible to transmit them successfully. 



/I(h.u) OF C/down 



Fig. 4: The data flow of a link. 


For downstream, the gateway is the only node that au¬ 
tonomously generates packets with rate gdown- The respective 
receivers v of a link (u, v) transmit traffic to be forwarded to 
their proper descendants with rate /i( u> „). This traffic is split 
up between the outgoing links. Each link ( v,w ) going out 
of v gets that fraction of the traffic Xr ViW \ that corresponds 
to the number of nodes reachable via this link (1 + 7™) in 
relation to the total number of nodes reachable via v, namely 
7 v . The second difference to upstream is the fact that each 


node will consume the fraction of the traffic intended for this 
node. Altogether, this results in the recursive expressions 


l + 7-u; 
7t> 

l+7u 


1 ffdo 


v = r A (v, w) G L 


down 


7v M(u,v) T (u,v) A (v,w) G ■^down 5 (17) 

9up + Yl(u, v ) e L up B(u,v) («iw) £ L up 


/-t( V,W ) ’ P(v,w) 


l+7w 

1 


(tl, tit) G Ldown 
(v, w) G L up 


• (18) 


Since the probability of a pending packet is modeled as 
Poisson distributed, it is given by qi = 1 — e~ Xl . 


V. Model for IEEE 802.15.4 CSMA/CA 

The IEEE 802.15.4 wireless standard uses CSMA/CA for 
media access control. There are other options, such as a beacon 
enabled mode and guaranteed time slots (GTS), but they are 
not easily realizable in multi-hop networks. 

The MAC layer is modeled by a Markov chain that is 
depicted in Fig. 5 with the following states. 

Sidle No transmission pending. 

Si,k,j In the i th backoff stage after j transmission attempts. 
CCA in k steps. 

0-i,qj During a successful transmission after q time steps 
and j preceding transmission attempts. 

0-2, q ,j During a colliding transmission after q time steps 
and j preceding transmission attempts. 

P(£ —» Q is the probability of going from £ to ( in one step. 
The probability that the sender will leave the idle state within 
one time unit is given by qi as calculated in the previous 
section. Therefore, the probability of staying in the idle state 
is 


P(0 id i e -> Oidle) = 1 - qi- (19) 

Otherwise, the MAC layer starts the backoff period. That 
is it waits for an integer random time span, m is the con¬ 
stant parameter macMaxCSMABackojfs that corresponds to 
the maximum number of backoffs as defined in [6], Note 
that, according to the standard, m does not include the first 
unconditional backoff, so the actual number of backoffs is 
to+ 1. The maximal backoff time span W, — 1 increases with 
the number i G [0, m] of failed channel accesses, with 

( 2 m ° i = 0 

Wj = < 2* • Wo 0 < i < m = mb — mo , (20) 

{2 m -W 0 i>m 

with the constant parameters too and to b that are the initial 
and maximum backoff exponents. The die is symbolic for the 
associated choice. It is not a state in itself, but just bundles the 
different ways to go into the backoff period. The probability 
of choosing a specific state is W t , therefore 

P(0idle^0o,k,o) = ^ r , 0 < fe < Wo — 1. (21) 

Wo 
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Fig. 5: Sketch of the IEEE 802.15.4 CSMA/CA technique. 


Within one time unit, the backoff duration is decreased by one 

P(9i,k+l,j “A = lj (22) 

0<i<mAO<fc<Wj — 1 A 0 < j < n. 


As soon as the state 9i : o,j is reached, channel sensing takes 
place. If the channel is sensed busy with probability 1 — cq as 
calculated in Sect. VII, the channel access retry counter i is 
incremented and the backoff procedure is started again 


P(0i-l,O,j ~A 0i,k,j ) — TTT 


(23) 


0<i<mA0<k< Wq — 1 A 0 < j < n. 


If i has reached m, the packet will be dropped 


P(9m,0,j ^ 9i d [ e ) — Ct/: 


0 < j < n. (24) 


Otherwise, if the channel is sensed free, the transmission 
takes place. If the packet or the associated acknowledgement 
sent out by the receiver will eventually collide, 0_2,o,j is 
chosen as next state 

P(9i,0,j “A 0-2, 0,j) = (1 — a l ) ' PnoACK,Zi (25) 

0 < * < to A 0 < j < n. 

Of course, the real hardware does not know if the transmission 
will collide, but in the analytical model the probability P no ACK,; 
can be computed from the collision graph as presented in the 
next section. In case of a upcoming successful transmission, 
the 0_i,o,j state is chosen 

P(0j,O,i -A 0-1,0,i) = (1 — Oil) ■ (1 — PnoACK, ;), (26) 

0 < i < m A 0 < j < n. 


During the transmission, the counter h is increased until the 
time is over 


P(9-i, h j -A 0_ i, h +i,j) = 1, 0 < h < L c A 0 < j < n, (27) 
P(0— 2 ,/ij 9-2,h+i,j) = 1,0 < h < L s A 0 < j < n. (28) 

After a collision, the retransmission counter j £ [0, n], with n 
being the maximum number of retries macMaxFrameRetries , 
is incremented, the channel access retry counter i is reset to 
zero and a new backoff period starts 

P(0— 2 ,L c ,i— 1 -A 00 ,k,j) = 0 < j <n. (29) 

w 0 

Finally, the idle state is reached after a successful transmission 
or as soon as the maximum number of retransmissions is 
reached and the packet is dropped 

P(0-l,i s ,j —> 0 idle) = 1, 0 < j < 71, (30) 

P(0— 2,L c ,n -A 9 ld le) = 1. (31) 


All transition probabilites not specified before are zero. The 
stationary distribution of the Markov chain yields the prob¬ 
ability that a given link is in a specific state as derived in 
Appendix A. In particular, the probability of being in a generic 
sensing state is given by 


n 


1 - n m+1 
_u 1 a l 

— 00,0,0—j- 

1 — Oil 


1 - Vi +1 

i -yi 


(32) 


with the shortcut y t := P n0 ACK,i • (l - ol™ +1 ) and 6 0 ,o,o, the 
probability of being in the first channel sensing state, given by 


-i i i /o \min(m,m)+l q min(m,m)+l 

‘ =1 UiAEiA- + ; 

oo,o,o 1 — 2 • ai 1 — on 


(2 mb + l)a™ +1 ^1 — a 


max(0 ,m—my 


l - 


Vi 


n+1 


1- ai 


1-2/; 


+ a-«r +i ) 


m+l^ 1 2 h 


,n+1 


1 - 2 /; 


-(Lb ’ (1 — PnoACK,;) + L c ■ P no ACK,;) 



1 - y? +1 ( m+1 , 1 - PnoACK,; \ 

!-2/; ^ 1 (l-ar +1 )'V 


• (33) 


VI. Simultaneous Transmissions 

The probability that the sender of a link j £ L does not 
start a transmission at a given point in time is given by 

TjOlj + (1 ~Tj). (34) 

with tj, the probability that the sending node of j is trying 
to access the channel and ctj, the probability that the channel 
would be sensed busy (see Sect. V). So either the sender tries 
to access the channel, but senses it busy ( Tjctj ) or it does not 
try to access the channel at all (1 — Tj). For a given set of 
links S C L, the probability that at least one j £ S starts a 
transmission is therefore given by the complementary event of 
the event that all nodes are not sending 

Q(l> S) = 1 — P[ (TjOtj + (1- — Tj)). (35) 

je s 
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The same expression can be derived from the expression 
Pr[_4;] in [2] as conducted in Appendix B, but is less complex 






















and takes much less computation time. For time intervals of 
t > 1, we take into account that the random variables are 
Poisson distributed, so 

Q(l, S) = 1 — e~ a ^ e~ a = 1 - Q(l, S) (36) 

for some rate er. For larger intervals, the rate increases pro¬ 
portional, so the probability of at least one transmission start 
within this time interval is 

Q(t, S) = 1 - e~ a t (37) 

= 1-(1-Q(1,S))‘. (38) 

VII. Analyzing the Neighborhood 

The formula derived in the previous section can be used to 
calculate collisions that occur if at least two transmissions ar¬ 
rive at the same receiver at the same time. In order to represent 
packets (sent along a link) as well as acknowledgements (sent 
in the reverse direction), even under the hidden node problem, 
an enhanced conflict graph is introduced that considers all 
constallations of senders (in the following: sends a packet, 
but receives an acknowledgement) and receivers (receives a 
packet, but sends an acknowledgement). 




Fig. 6: The four basic relations for conflicts of two links. 

For a link (vi,Wi) G L the possibly disturbing links 
(v 2 , W'i) £ L are expressed in terms of the basic constellations 
as depicted in Fig. 6 and formulated as 

Ss,(v i,®o = {(^ 2 , w 2 ) G L | D(vi,v 2 ) A vi ± v 2 }, (39) 

that is the sender of one link receives transmissions from the 
sender of another link, 

= {(v 2 , w 2 ) G L | D(wi, v 2 ) A vi ± v 2 }, (40) 

that is the receiver of one link receives transmissions from the 
sender of another link, 

Stz,(vi, W i ) = {{V2, w 2 ) G L | D(V!,W 2 ) AV!^ v 2 }, (41) 

that is the sender of one link receives transmissions from the 
receiver of another link, 

'R-n,(v i,»i) = {{V 2 ,w 2 ) G L | D(w 1 ,w 2 ) A vi ± i> 2 }, (42) 

that is the receiver of one link receives transmissions from the 
receiver of another link. 

In the following, CPi denotes the event that a packet 
collides on a link l, while CAi denotes the event that an 
acknowledgement collides. The constellation that the sender 
v 2 of the link j = (v 2 , w 2 ) is in the range of both the sender 


v\ and the receiver w 1 of a link l = (vi,W\) can be expressed 
as 

j g Hs,i r\Ss,i. (43) 

In this constellation, both senders will wait for each other if 
they sense an ongoing transmission. Although, a sender might 
sense within the turnaround time of the other one and vice 
versa and they therefore start transmitting within an interval 
of 2 time units. So the probability of a collision of two packets 
on links that obey this constellation is expressed as 

P{CPi,o) = Q(2,n s ,inSs,i). (44) 

If receiver Wi might receive a packet from sender v 2 , but the 
transmission of v 2 is not recognized by vi, the channel sensing 
has basically no effect. This is called the hidden node problem. 
Two transmissions of length L p might overlap in an interval 
of 2 L p . This is calculated as 

P(CPi, i) - Q(2L P , H s ,i \ S s ,i). (45) 

The transmission can not only collide with another packet, 
but also with an acknowledgment. If all involved nodes are in 
range, the new packet might collide with an acknowledgement 
if the carrier sensing takes place between the arrival of the 
original packet and the beginning of the acknowledgment 

P(CPi, 2 ) = Q(i,s Stl ns n>l mZn tl ). (46) 

If the senders can not hear each other, the sensing might also 
take place at the end of the packet, 

P(CPi, 3 ) = Q(2, (S n ,i n TZ n ,i) \ S s ,i). (47) 

If it were to occur earlier, the transmission would collide 
with the packet, so it should not be counted here, too. If the 
acknowledgement can not be heared, the packet transmission 
might take place during the whole acknowledgement 

P(CPi A ) = Q(L ACK , (S S:l n K n ,i) \ S n ,i). (48) 

This case is extended in the following constellation by the 
sensing before the pause (cf. P(CP 1 , 3 )) 

P(CP 1 , 5 ) = Q(L ack + 1, (Us, I n Hu, 1 ) \ Ss,i \ Sk,i). (49) 

Finally, the packet might overlap with an acknowledgement 
during the whole transmission, if the sensing of both senders 
is not effective 

P(CP 1 , 6 ) — Q(L p + L ack ,1Zk,i \ Ss,i \ $k,i \ (50) 

In addition, the acknowledgement itself might be affected by 
collision with a packet 

P(CAi fi ) = Q(i,Ss,inn s ,i), (51) 

and in particular if the acknowledgement can not be heared 

P(CAi A ) = Q(L ack , S sj \ H s , t ). (52) 

The probability that a packet collides with at least one other 
packet or acknowledgement is 

P(CP l ) = p( [ \JCP l }j. (53) 

Note that the events are not mutually exclusive, so inclusion- 
exclusion principle has to be applied for the calculation. 
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Furthermore, a transmission of b bytes might be dropped 
because of a bad link with probability PER;, / as given in (4). 
Taking this into account, the overall probability that a packet 
is not successfully received is 

P(CPi) = P(CPi) + (1 - P(CPi)) ■ PER b ,/. (54) 
The combined probability of a lost acknowledgement is 

P{CAi) = p[{jCA l )j. (55) 

Taking a possible bad link into account, the probability of 
losing the acknowledgment is 

P(CAi) = P(CAi) + (1 - P(CAi)) ■ PER Back , ; . (56) 

Finally, the overall probability that no acknowledgement ar¬ 
rives at the sender is 

PnoACK,* = P(CPl) + (1 - P(CPl)) ■ P(CAl). (57) 

The probability ai that the channel is sensed busy is 
calculated along the same line as the probability that either 
a packet transmission or an acknowledgment arrives at the 
sender 

&pkt,i — ,i) i (58) 

&ackf ^ 2 (-^ ack > &lZ,l)i ( 59 ) 

&pkt,l T &ack,l &pkt.,l * (%ack,L• ( 60 ) 

Note that all sets can be computed offline, so the actual 
computation is linear with the number of neighbours. 

VIII. Link Reliability Considering 
Retransmissions 

It turned out that the reliability of packet transmission is 
not independent of the current retransmission attempt and this 
effect has a major impact on the results [1], Therefore, the 
model contains an elaborate handling of this matter. 

The probability of a mutual disturbance of at least two 
hidden nodes is calculated 

P(CB la ) = Q(2L p + 2, {Us,i n S n ,i) \ S s ,i). (61) 

After such a mutual disturbance took place, those nodes will 
issue a retransmission. The backoff exponent is reset to the 
initial backoff counter mo, so there are Wq = 2 m ° possible 
backoff time spans. The event that the retransmission collides 
again, given i other nodes with pending retransmissions is 
denoted as CIZi , 2 ■ For two nodes the probability is 

u> = min(Wo — L p — 1,0) 

i\ i ' 2 h i ui + w 2 (62) 

= w o 

Of course, this expression is only a lower bound of the 
actual probability that there will be a repeated collision. For 
example, there might be several initial collisions at the same 
time leading to a higher probability of a repeated collision. It 
can be calculated by summing over all possible combinations 
of mutual disturbance. However, the price for this gain of 
accuracy is too high considering the exponential growth in 



computational complexity, in particular for large scale net¬ 
works. Analogous to (44), even senders which can mutually 
sense their transmissions might be affected by packet collision. 
The probability that at least two senders which can mutually 
sense their transmissions are affected by mutual disturbance is 

P(CB lt i) = <2(2, 71s,! n S KJ n Ss,i). (63) 

CJZi,i J is defined analogous to CR,/,-/ and the corresponding 
probability of a repeated collision of two nodes is calculated 
as 

P(CJh/) = WV 1 . (64) 

With these quantities, an absorbing Markov chain is built 
with the states 

k succ Successful transmission 
K c f Channel access failure 
k o,o No preceding mutual disturbance 

o Hidden node(s) with pending retransmission 
/^o,l Visible node(s) with pending retransmission 
Ki,i Hidden node(s) and visible node(s) with pending 
retransmission 

In the following, all unspecified transition probabilities as well 
as P(CRi, 2 °) and P(CRi, i°) are zero. k succ and n c f are the 
absorbing final states, so 

P^^succ ^ tv succ ) 1? (65) 

P{Kcf Kcf) = 1. (66) 

The probability of a channel access failure is always the same 
P(k p ,,-> K cf ) = a™ +1 . (67) 

All following transitions need to take this into account with 

a = i -«r +1 - (68) 

The probability of a successful transmission is given by the 
probability that neither a conventional collision takes place 
with P(jCPi) nor a repeated collision 

P(K p , q —> K succ ) = A ■ (1 — P{£Pl U CRl,2 P U CRl,i q )). 

(69) 

Note that P(CPi ) is used instead of P n 0 ack,z< because even if 
the acknowledgement does not arrive at the sender, the packet 
itself might be transmitted successfully. The transmission to 
Ko,o takes place if no repeated collision takes place, but the 
packet collides anyway with another transmission 

P(Kp, q —f K 0 ,o) = A • (P(CPi) - P(C-B lp _ U CB U i)) (70) 

• (l-P^/UCRu’)). 

All remaining state transitions indicate that any repeated col¬ 
lision takes place 

P(Kp, q —f Ki,o) = A • P{CB t , 2 U Clh,2 P ) (71) 

• (1 — P(CBi,i Uffi,/)), 

P(Kp, q —f Ko,i) = Pi ■ (1 - P(CB la u CRi, 2 P )) (72) 

■P(CB L1 U CR ul q ), 

P(K P , q —f ki,i) = /3i ■ P(CBi, 2 U CRi/) (73) 

• P(CB,i U CR^). 

The probability Ri that a packet was successfully transmit¬ 
ted is described by the probability of reaching k succ from ko,o 
in at most n + 1 steps, calculated by the corresponding power 
of the transition matrix. 
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IX. Path Reliability 


The end-to-end reliability in upstream direction lR uPi „ is 
the probability that a packet sent by a node n arrives at the 
gateway. It is given by the product of the reliabilities along 
the routing path 




up,c — 


f 9^up,p * R ( c,p) c 7^ t A (c, p ) £ L up 
\l c = r 


(74) 


By construction, there is always exactly one p with (c. p) £ 
L up unless c is the gateway. The reliability in downstream 
direction (Kdown.n* that is the probability that a packet sent by 
the gateway arrives at its destination, is given correspondingly 
by multiplication of the reliabilities in downstream direction 



f -Tidown,p ’ R(p,c) C 7^ T A (p, c) £ Ldown 
\l c = r 


(75) 


List of Symbols 

ai Probability of sensing a busy channel. 

&pkt,i Probability of sensing a packet. 

cx a ck,i Probability of sensing an acknowledgement. 

B Bytes per packet. 

b^ Probability of being in state 9^. 

Back Bytes per acknowledgement. 

BLR/ Bit error rate. 

/?; Shortcut for 1 — ai m+1 . 

CAi Event that an acknowledgement on l collides. 

CBi Event that transmissions mutually collide. 

CPi Event that a packet on l collides. 

CJZi Event that transmissions collide again. 

D(v,w) Predicate that describes potential collisions be¬ 
tween nodes. 

7 n Number of proper descendants of node n. 

Sdown Packet generation rate in downstream. 
p up Packet generation rate in upstream. 

Tdown Packet interval in downstream. 

IFS Inter-frame spacing as defined in the standard. 

/ Llp Packet interval in upstream. 

k s State in the retransmission Markov chain. 

L p Duration of packet sending. 

L The set of active links. 

L A c K Duration of acknowledgement. 

LAi Event that an acknowledgement on l gets lost. 

L c Total time for an unsuccessful transmission. 

Ldown The set of active links in downstream direction. 

A/ Packet sending rate. 

CPi Event that a packet on l gets lost. 

L s Total time for a successful transmission. 

L up The set of active links in upstream direction. 

to Maximum number of channel sensing attempts 

( macMaxCSMABackoffs ). 

Too Initial backoff exponent. 

mb Maximum backoff exponent. 

to Shortcut for nib — too- 

pi Packet forwarding rate. 


N Total number of nodes in the network. 

n Maximum number of retransmission attempts 

(mac Max Pm me Retries). 

Pdist Minimum power that interfers with a reception. 

PER/, / Packet error rate for a transmission of b bytes. 

P N Noise power. 

Puoack.z Probability that no acknowledgement arrives after 
sending a packet over l. 

P rx Reception power. 

P tx Transmission power. 

qi Packet sending probability. 

Q(t. S) Probability that at least one link in S starts a 

transmission within t time units. 

r The gateway node. 

(Ttdown,; Reliability that a packet sent by the gateway 
arrives at n. 

Ri Reliability of packet transmission on a link. 

Tln.i Set of links with receivers in the interference range 
of the receiver. 

1Zs,i Set of links with senders in the interference range 
of the receiver. 

91 U p,z Reliability that a packet sent by n arrives at the 
gateway. 

Si, Backoff unit time (used as basic time unit). 

SNR; Signal to noise ratio. 

S-jzj Set of links with receivers in the interference range 
of the sender. 

Sg.i Set of links with senders in the interference range 

of the sender. 

tack Time interval between packet transmission and 

acknowledgement. 

r; Probability of a channel sensing attempt. 

tm,ack Time to wait for a lost acknowledgement 
( macAckWaitDuration ). 

9 S State in the MAC layer Markov chain. 

T(p, c) p is parent of c in the routing tree. 

Wi Maximum backoff time span for the ; lh channel 

sensing. 

yi Shortcut for P noA CK,i ■ (l - ct; m+1 ). 
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Appendix A 

Stationary Distribution 


In the asymptotic case, the probability b^ of being in state is calculated by summing over the probabilities of the preceding 
states S, weighted by the transition probability. 

&C = • P ( 0 « “►*<) (76) 


Therefore, from (21) and (22) results 

bo,k ,0 = 1 -jrbidle + fro,fc+l,0> 

IVo 

Since the maximum backoff period is Wi — 1, it holds 

h,Wi,j = 0 

so the recursion resolves to 

fro,0,0 = Wo - -jp bidle = qi ■ bidle- 
W 0 


The same for (23) results in 


on 


frz,0,j fri—l,0,j — tit ‘ fri—1,0,7? 

= Ot] ■ fro.O.j 


and for (29) 


fro.O.j = Wo ' ^rfr-2,L c ,j-l = b-2,L c ,j-l, 


together with (25) 


fro.o ,j = ^ fri,0,i-l • (1 — Oil) ■ P no ACK,Z 
i=0 

m 

= bo,0,j-l ’ (1 — Oil) ■ PnoACK.Z ' Ctj 

i=0 

/ m \ 

(1 — Ctj) • P no ACK,Z • ^ «£ 


= fr( 


' 0 , 0,0 


i=0 


and taking the sum as the partial sum of the geometric series 

/ 1 _ \ 

= fro,0,0 • (! — «;)• PnoACK,i 


1 - cxi 

= fro,0,0 • (PnoACK.Z • (l — a ; m + 1 )) J 
= fro, 0,0 • 3/f - 


0 < fc < W - 1. 


0<*<mAO< j<n, 


0 <i <m— 1 A 0 < j <n, 


0 < j < n. 


For each k G [0, W - 1] A * > 0 holds with (21), (22), (23) and (29) 

for i = 0 A j = 0 

bi,k,j — fri,fc+lj 


w 0 Uldle 

b ■ , 
w *-l 


for 0 < i < to — 1 A 0 < j < n 
^b- 2 ,L c ,j-i for i = 0 A0 < j < n 


and taking (79), (80) and (82) into account 


Qi 1 u 
Wo qi W°’0 
ai 1 1-. 

Wi a t 

Wo fe °-°d 


for i = 0 A j = 0 

for 0 < * < to — 1 A 0 < j < n 

for * = 0 A 0 < j < n 


— bi,k+i,j + bi,o,j ■ 


(77) 

(78) 

(79) 

(80) 
(81) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 


(90) 

(91) 



The recursion dissolves to 


bi.k.i — b' 


i,k,j — ^2,0 ,j 


Wi - k 
W ’ 


and together with (81) and (88) 


bi,k,j = b 0 , o,o -Vi-ot]- 


W z -k 

m ■ 


The probabilities in the stationary distribution have to sum up to one (normalization condition), that is 


m n Wi — 1 


n /L s —1 


2—0 j— 0 k— 0 


EE E bi,k,j + EE b-l,h,j + E] b-2,h,j I + bidle — 1- 


j =0 \ h =0 


From (93) results 

m n Wi — 1 


m n Wi — 1 


E E E b kkj - E E E 


o • y\- a\ ■ 


Wi - k 


2—0 j—O A:—0 


2—0 j— 0 A;—0 


— &o,o,o EE^°i 

2—0 j =0 
n m 

= bo,o,o E yi • E a ? 




Wi + 1 


J—0 2=0 


splitting the second sum at m 


n / min(m,m) 


Wi + l 


W+ +1 


i 2 i Wo + 1 


i 2 m W / o + 1 


= Vo,oE^' E + E 

j=o y 2=0 2=m+i y 

„ /minhn.m) Q + 1 - 2*W 0 + l\ 

= bo,o,oJ2y l \ E a r — 2 —+ E - 2 - 

j—0 \ 2=0 2=m+l y 

, 1 _ n+1 / min (m,m) m 

- ° 2 °’° E aj-(2W 0 + l)+ E a* • (2^2 m ° + l) 


h 1 „"+! / min( 

On.o.o ± ~ Vi / v 


fro,0,0 1 ~ Vi 
2 1 - Vi 


E (2a;)* • Wq + E a) + (2 mfc - mo+mo + 1) E 


&o,o,o 1 - vr 
2 1 2/; 


,n+l /min( 


E (2a;)* • Wq + E a) + (2^ + 1) E 


m+1 i—m— 1 


; 1 „ n+1 / 1 CO \min(m,m) + l 

&o,o,ol~y; |^. l-(2aO 


2 1 - 1 /; 


1 - 2a; 


min(m,m)+l 

1 - a; 


+ (2 mb + l)a)” +1 E 


&o,o,o 1 - V? +1 ( u/ 1 - (2a;) min(m ’ m)+1 ^ 1 - ar (m ’ m+1) 

° ^2a; 


(2 m! > + l)a™ +1 (l 


1 - a; 


It follows from (26) and (27) 


b-l,h,j = b-l,o,j = E(! _ ‘ (1 — -PnoACK,;) • bi, 0 


— (1 — -PnoACK,;) • (1 — a;) • 6o,0,0 ■ y{ • E 

i=0 

= (1 - PnoACK,;) • bo,0,0 ■ yj • (l — a T +1 ) 
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and analog from (25) and (28) 


b-2 ,h,j = PnoACK,/ ’ bo,0,0 ’ Ul ' (l _ a T + ± ) 


— yj +1 • 6o,o,o- 

Inserting these into the second term of the normalization condition gives 

n /L s — 1 L c —1 \ n 

y (y b-i th j + y b- 2 ,h,j ) = y( L ° • (o _ aoack,o • (>0,0,0 • $ ■ (1 - 


j =0 \ h =0 


+L C ■ ^PnoACK, Z • (>0,0,0 ' yj • (l — a " , + 1 ))) 

n 

= bo, 0,0 • (l — a T +1 ) ' (L s ■ (1 — PnoACK,>) + L c ■ P no ACK,z) ’ ^ Ul 

3=0 
1 — 

= bo, 0,0 ' (l — «r +1 ) ' (-^s ' (1 — PnoACK,>) + Lc ' PnoACK,Z) ’ ~7- ~ 

1 - Vl 


The idle state probability is according to (19), (24), (30) and (31) given as 


bidle — (1 9z) ' bidle P ^ ' (o?Z * b rn o,j P P b—2,L c ,n 

3=0 

n 

bidie — (1 — Qi ) • bidie = ' ap • j// • (>o,o,o + (1 — PnoACK,;) • (>o,o,o • yf ■ (l — ar +1 )) P 


r +1 ))+y? +1 -bo, 0,0 


bidle = bo, 0,0 ■ y y 3 i («r +1 P (1 - PnoACK,0(l ~ C^ 1 )) P ^ 




(>0,0,0 (l_yy( a rn+l p (1 _ P BoACK l )(l - «p+l)) P y ; n+1 ) . 


Finally, inserting (104), (112) and (116) into the normalization condition (94) gives 


,-i ll-2/r +1 / TT , l-(2a J ) min( "’ 

bo,o,o — 2 l _ yi \ W ° 1 _ 2ai 


min(m,m) + l min (m,m- 

, 1 ~a l 


+i) ( 2 mt + l)a” 1+1 (l 


max(0 ,m—m) 


1 - ai 


1 - ai 


1 — v n+1 

P (l — < + 1 ) ' (P« ‘ (1 — PnoACK,>) P Pc ' PnoACK,z) ' ~J^- 

+ \i ' ( y " +1 + 1 ~ ryr ^ +1 + (1 - p "°ack,z)(i - «r +1 ))) ■ 

The probability 77 is then calculated by summing over all bi o,j with 0 < i < m and 0 < i < n, giving 

m n 

t i= yy ko ,3 

2=0 j—0 
m n 

=EE a\bo,o,j 

2=0 j = 0 

m n 

= (> 0 , 0,0 E yf 

2=0 J =0 

1 - «r +11 - %” +1 

— 0 0 , 0 , 0 —;- 7 —• 

1 - ai l-yi 
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Appendix B 

Alternative Derivation of <3(1, S) 

The expression in [2] iterates over all combinations of one up to all neighbouring links l £ S. This can also described 
power set of S excluding the empty set 

P(S) \ 0 = {M\M CSAM^I}. 

With this, the expression can be written as 

qu,s)= ^ n a-Jfi-n*) 

MeV( S)\0 ViGM / \ieS\M ) V ieM ) 


= e iH n (i-^)i(i-n 

M6P(S) ViGM / \ies\M 


n r - nt 1 -.) En- 


a.i I - I — 77) I I l — II a,; 

ieM ) ViG0 / ViGS / \ iG0 


= e n r * n a-^Ki-n ^)- 0 

JU6P(S) \i£M / VieS\M 


iGM 


= e f n r * n n ) - e f nH n - t *) ) n 




Mev(s) \ieM / Ues\M 


MG'P(S) \iGM / \ieS\M 


ieM 


e (lm n - e (n 

Mev( s) ViGM / \ies\M 


T i a i II n a-*o 

Mev(s) \ieM ) \ies\M 


Note that for j £ S holds 


P(S) = {M U {j} | M £ V(S \ {j})} U {M | M £ P (S \ {j})}. 

Therefore, the second term can be split up 

e (n r ^) (n ^ 

MGP(S) ViGM / \ieS\M 

= e ( n ( n - < 


- TqCMq 


e n ^ 1 n 


— T i) 


MeV(S\{j}) \.GM 


oe(s\{f})\M 


+ ^ n ^ n (! - t k 

MgP(S\{j}) ViGM / \ieS\M 

(! - Tn) yi (n nai ) ( n - r *) 

MgP(S\{j}) ViGM / V*G(S\{j})\M 


/ 


= ( r 3 a 3 +(! - ^)) e n Ti<xi ) n (! - T i) • 

MeV(S\{j})\ieM J \iG(S\{j})\M / 

This step can be recursively repeated for all nodes in S, finally resulting in 


E fn 


Tiai ii n ( x - r E= n < ' T :i a i +( x - t j))- 

MeV{S) \i£M J \ieS\M 


l GS 


The same calculation for the first term in (127) results in 


Finally, the result is 


e iH n ( i-t *) ) = n( T J + ( i - 7 j)) =i - 

MeV( S) ViGM / ViGS\M / jGS 


Q(l, S) = 1 - YiiTjCXj + (1 - Tj)). 
fes 


as the 

( 122 ) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 

(129) 

(130) 

(131) 

(132) 

(133) 

(134) 

(135) 
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